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Abstract. In bipartite quantum systems commutation relations between the 
Hamiltonian of each subsystem and the interaction impose fundamental constraints 
on the dynamics of each partition. Here we investigate work, heat and entropy 
production in bipartite systems characterized by particular commutators between 
their local Hamiltonians and the interaction operator. We consider the formalism of 
[Weimer, EPL, 83:30008, 2008], in which heat (work) is identified with energy changes 
that (do not) alter the local von Neumann entropy, as observed in an effective local 
measurement basis. We demonstrate the consequences of the commutation relations 
on the work and heat fluxes into each partition, and extend the formalism to open 
quantum systems where one, or both, partitions are subject to a Markovian thermal 
bath. We also discuss the relation between heat and entropy in bipartite quantum 
systems out of thermal equilibrium, and reconcile the aforementioned approach with 
the second law of thermodynamics. 


1. Introduction 

The emergence of thermodynamic behaviour within quantum mechanical systems has 
attracted much attention in recent years This interest is motivated by the desire 

to gain a deeper understanding of thermodynamics, as well as the need to obtain 
generalizations of thermodynamic concepts and constraints to quantum settings. A 
conceptual understanding of work, heat and the second law in the quantum domain has 
therefore been central to research in quantum thermodynamics PHEIEIITIIH]. Owing to 
advancements in quantum technologies, experimental tests of these concepts are within 
tangible reach p HQl HH [12]. The developments in this area allow investigation of the 
relevance of coherence in the efficiency of quantum engines PI and can guide the design 
of nano-mechanical devices with thermodynamic functionalities danni. 

The conventional dehnition of work and heat for a quantum system evolving under a 
time-dependent Hamiltonian considers the change in the internal energy (7 of a system as 
U = Tr{pi7}-|-Tr{pi7} and identihes the hrst (second) term as the work (heat) flux |16j . 
This division assumes that a classical, external, driving gives rise to a time-dependent 
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Hamiltonian. In contrast, we are interested in autonomous bipartite quantum systems 
in which the inherent quantum mechanical interaction between the two partitions results 
in internal transfer of energy. This internal transfer can subsequently be identihed as 
either a heat, or a work, flux [m ng [IS]. There are plenty of quantum settings where 
such is the case, but an interesting example is a prototype biological dimer where the 
interaction between electronic and intramolecular vibrational motions results in non- 
classical behaviour and internal energy transfer [20] . A framework for dehning work and 
heat fluxes in autonomous quantum systems was presented by Weimer and colleagues 
in Ref. [T7]. This approach, henceforth referred to as the Weimer framework or the 
energy flux formalism, argues that time-dependent Hamiltonians arise from tracing out 
the degrees of freedom of the driver and need not have a classical, or external, origin. 
By dehning an effective local energy basis, heat (work) hux is identihed with energy 
changes that (do not) alter the the local von Neumann entropy. This framework is of 
considerable conceptual interest as it aims to extend the connection between entropy 
and heat to hnite, out of equilibrium quantum systems. Despite its conceptual ease and 
formal structure, this approach has not been widely used in the investigation of work 
and heat in quantum systems EH- This is arguably due to the diherences between the 
work hux dehned in this way and the extractable (measurable) work; the latter being 
the focus of most prior research mu- As work done on a system is not a quantum 
observable [5], its quantum statistical properties and the backaction of the measurement 
apparatus need to be considered within any formalism. These issues have partially been 
addressed within the energy hux formalism ra but remain under scrutiny. 

The simplicity of the energy hux formalism, however, merits the investigation of 
aspects scarcely explored such as the implications of conservation laws on the nature of 
energy exchanges and how the framework can be generalized to open quantum system 
dynamics. Conserved quantities (dehned by Hermitian operators that commute with 
the total Hamiltonian) and their associated symmetries are of broad interest in physics. 
Moreover, transformation of symmetry-imposed constraints in open quantum systems 
has gained renewed attention in recent years [22]. In this article, we discuss work 
and heat exchanges within bipartite systems characterized by particular commutation 
relations between the local Hamiltonians and the interaction operator. We show 
how conserved quantities associated to a given partition, that arise from particular 
commutation relations, give rise to an asymmetric structure for the work and heat huxes 
of the subsystems. We furthermore demonstrate the consistency of the approach when 
one partition acquires the character of a macroscopic bath, and extend the formalism 
to situations where one, or both, partitions are subject to a Markovian thermal bath. 
The connection between heat, as dehned in this framework, and entropy is explored in 
the last section. This connection allows us to discuss the second law of thermodynamics 
within the energy hux formalism and demonstrate the compatibility of the approach 
with the second law. The present work strengthens the energy hux formalism and 
provides further evidence that thermodynamic quantities such as work and heat can be 
generalized to quantum systems far from thermal equilibrium. 
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Figure 1. A schematic illustration of the bipartite quantum system investigated in this 
article. In section 2, we focus on the closed bipartite system (without a heat bath). In 
section 3, partition I is subject to a thermal bath. In section 4, both partitions interact 
with a thermal bath. In our notation the total heat flux into system I is written as 
Qi = Q 2 -^l + Qbath A-i-i where Q 2 -^l and Qbath A-i-i are the heat fluxes from system 
II and the bath A respectively. The work flux into system I is identified as Wi. 


2. Work and heat in closed bipartite quantnm systems 

2.1. Review of the concepts 

In this section we briefly review the division of the energy exchange between bipartitions 
of a quantum system into heat and work, as put forward by Ref. na. Within the 
framework open quantum systems, the total system is denoted as closed when it does 
not exchange energy or matter with an environment and is described by the Hamiltonian 
H = Hi + H 2 + Hi2 . Here Hi and H 2 are the individual Hamiltonians of subsystems 
I and H, and H 12 specifies the interaction between the two partitions. Notice that in 
classical thermodynamics our total system would be referred as isolated. The density 
matrix of the total system can be written as 

p{t) = pi{t)^P2{t) + Ci2{t), (1) 

where pi{t) = Tr 2 {p(t)}, P 2 {t) = Triand Ci 2 {t) represents the correlations 
between the two partitions. We assume an initial separable state of the form p(0) = 
Pi(0)<8)p2(0), with no bipartite correlations. Starting from the equation of motion for the 
total system, and taking the trace with respect to system H, the Liouville-von Neumann 
equation for system I can be written as 

pi{t) = -i[Hi + Hf{t),pi{t)]-ii:i2{[Hi2^Ci2{t)]], (2) 

where the effective Hamiltonian of system I is identified to be 

Hf{t) = Ti2{Hi2[li^P2{t)]}. (3) 

The energy changes need to be defined with respect to a local effective measurement 
basis (LEMBAS), which is determined by the experimental apparatus. This formalism 
assumes a LEMBAS that includes both the eigenstates of the bare system Hamiltonian 
Hi, and the shift in the eigenenergies of system I due to the interaction with system 
H. Since Hi and Hi{t) may not be compatible observables at all times, the effective 
Hamiltonian is divided into two contributions 

Hf(t) = i/if (i) + i/if (0, 


( 4 ) 
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where (t) and are components of the effective Hamiltonian that commnte 

and do not commnte with Hi respectively. The commnting contribntion of the 
interaction induces a shift in the local energy of system 1. The local basis is thus 
associated to measurements of the following operator 

H[(t) ^ Hi + (5) 

and the internal energy change of system 1 can be identified as 

dUi = TTi{pidH[{t)} + TTi{H[{t)dpi}. (6) 

We now divide the internal energy flux of system 1 into two contributions: the energy 
change that has no effect on the local von Neumann entropy, is labelled as the work flux 
Wi{t), whereas the contribution that alters this entropy is identihed as the heat flux 
Qi{t). This constitutes a statement of the hrst law of thermodynamics for partition I. 
The final expressions read 

»',(() = ( 7 ) 

Q.(«)=lV,{i/:(«)£'ff|p(()]}. (8) 

where = —iTT 2 {[Hi 2 ,Ci 2 {t)]}. Note that at this point there is no assertion 

that work should be something more useful that heat [21] . Importantly, entropy changes 
can arise without an associated heat flux as is discussed in section 4. Heat flux can, 
alternatively, be expressed as 

Qi{t) = (g) l2),Hi2]Ci2{t)}, (9) 

where the cyclic invariance of the trace has been used. Similar expressions can be written 
for system H if measurements are performed in the LEMBAS associated with H 2 . 

Equations (Ej) and (I9|) merit two important observations. Firstly, work flux now has 
an additional contribution due to the non-commuting part of the effective Hamiltonian 
in Eq. (]4|). This is in contrast to the conventional dehnition that identihes work solely 
with the term Tri|pi(iiJ{(t)} in Eq. ([6]). Secondly, the heat flux into system I vanishes 
if at least one of the following is satisfied: (i) the dynamics does not induce correlations 
between the two subsystems or (ii) H'i{t) commutes with the interaction Hamiltonian. 
Thus, bipartite correlations are necessary, but insufficient, for the existence of a heat 
flux. We will discuss in section 4 the implications of this feature for the possible link 
between heat and entropy production in autonomous systems. 

From the above discussion, it is clear that commutation relations between the dif¬ 
ferent energy operators shape, in a non-trivial manner, the work and heat exchanges 
between the quantum partitions. In the following sections we consider three possible 
commutation relations between Hi, H 2 and H 12 , and investigate their consequences on 
work and heat fluxes. Without loss of generality, we consider an interaction of the form 
Hi 2 = Ai ^ B 2 , where Ai acts on system I and B 2 on system H. The effective Hamilto- 
nian of system I is thus given by Hi{t) = Ai{B 2 ). All results can readily be extended 
to the more general interaction H 12 = A^® Bi. 
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2.2. Fully commuting interaction 

A fully commuting Hamiltonian satisfies the relation: [Ai,ifi] = [B 2 ,Hq\ = 0. Since 
[Hf^{t), Hi] = 0, we identify Hi^a = Hf^{t), and = 0. Work on system I is thus 
given by 

vri(() = = {a,){B2} = 0, (10) 

where we have used (M) = —i{[M,H]) for M = li ^ B 2 . Similarly, one can show 
that W 2 = 0. From Eq. (|9]) it can be shown that the exchanged heat also vanishes, 
Qi(t) = Q 2 (t) = 0 . The energy exchange is thus zero, as expected for a commuting 
interaction. 

An example of such a system is provided by the spin-oscillator Hamiltonian with a 
dispersive interaction: H = + zzala -|- gazO^a, where is the population difference 

operator, ujq is the energy gap between the ground and excited state of the spin, {a, al} 
are the annihilation and creation operators of the oscillator, v is the oscillator frequency 
and g quantifies the magnitude of the spin-oscillator coupling. 

2.3. Partially commuting interaction 

We define a partiality commuting interaction via the commutator properties: [Hi, AA = 
0 and [H 2 ,B 2 \t^ 0 i.e. observable Ai is conserved. Since [Hf^{t), Hi] = 0, we identify 
Hf^{t) = Hf^{t), and Hf^ = 0. Work on system I is thus given by 

w,{t) = Tiqi:r';ff(()p,(()} = {/!,)(*). (11) 

Work done on system H can be computed in a similar manner. The effective 

rr nn 

Hamiltonian of system H is given by H^if) = {Ai)B 2 . Unlike the previous case, (t) 

nn nn nn 

has no commuting contributions with H 2 . That is, hffa ~ 0; ^ib ~ 


work done on system H is thus given by 

W2{t)= -lAl2{[H2,{Ai)B2]p2{t)], 

= - iTri,2{[li ® H 2 , Ai (g) B2]pi{t) (g) p2{t)], 

= -{Ai){i32). ( 12 ) 

The general work flux relation holding for any partially commuting interaction is thus 
Wi{t) + W2{t) = t). (13) 

The heat fluxes into systems I and H are computed from Eq. ([9]), leading to 

Qi(t) = 0, (14) 

Q2{t) ^ 0, (15) 


since [{H[{t) (g) l 2 ),hfi 2 ] = 0 while [(li ® H 2 {t)), Hu] 7 ^ 0. Note that control of the 
system Hamiltonian without introducing heat is reminiscent of classical driving. System 
H thus acts as a classical drive for system I, while the converse is not true. Moreover, 
for this commutation scenario, heat flux is exclusively from system I to system H, 
and the existence of correlations between the two partitions determines whether or not 
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system II receives heat. This observation is a direct consequence of conservation of the 
observable Ai, and does not violate the conservation of energy. 

Under what circumstances do the correlations, and therefore Q 2 {t), vanish? Note 
that the commutator [Hi,Ai] = 0 , ensures that Hi and Ai are simultaneously 
diagonalizable. Ai can thus be expanded in the energy eigenbasis of Hi. Denoting 
the eigenstates of Hi by {|*)}, Ai can be written as Ai = where {oj} are a 

set of coefficients. The interaction Hamiltonian is thus given by 

Hu = '^ai\i){i\ B 2 . (16) 

i 

From this Hamiltonian one can construct an evolution operator and determine the 
dynamics of the density matrix p{t) for an initially separable state p( 0 ) = pi( 0 ) ( 8 )/O 2 ( 0 ). 
If the initial state of system I is an incoherent mixture of eigenstates pi(0) = Yhi 
the dynamics are given by 








,i{H2+aiB2)t 


( 17 ) 


The total state only remains separable if all but one of the coefficients q, are zero, 
i.e.iwhen system I starts in an energy eigenstate of Hi. We conclude that the dynamics 
do not correlate the two subsystems, i.e. Q2{t)=0, if and only if system I is initialized 
in an eigenstate. 

An example of such a quantum setting is the displaced spin-oscillator Hamiltonain: 
H = ^az+i^a'^a+gazia^+a), where Hi = ^(Tz and H 2 = I'a^a. Denoting the eigenstates 
of the spin by {|i)} where |i) = {|e), l^f)}, for the initial state p( 0 ) = |i)(i| (D ^ 2 ( 6 ), 
the spin and the oscillator evolve independently and no heat is being transferred to 
the oscillator. However, if the spin is initially prepared in a mixed state of the form 
pi( 0 ) = c\g){g\ -t- (1 — c)|e)(e|, the dynamics correlate the two subsystems, generating 
a heat flux into the oscillator. To illustrate the latter point, let us assume that the 
oscillator is initially prepared in a coherent state, while the spin is prepared in a mixture 
of eigenstates. This example was considered in Ref. [I9] and the authors conclude that 
there is no build up of correlations between the spin and the oscillator. In the remaining 
part of this section, we demonstrate that this conclusion is incorrect; although the heat 
flux into the spin is zero (Qi = 0 ), there is hnite heat flux into the oscillator {Q 2 7 ^ 0 ) for 
an initial mixed state of the spin. The state of the spin and the oscillator may therefore 
become (classically) correlated. For this initial state, the work flux into the spin can be 
shown to be 


Wi{t)=g{l-2c){x{t)), (18) 

where x = a-|-ah We assume an initial coherent state denoted by |ao) with displacement 
xq = Re[ao] and zero initial momentum, po = Im[a!o] = 0 . The evolution of the density 
matrix of the oscillator is given by 

pit) = c|a_(()){a_(()| + (1 - c)|a+(())(a+(()|, 


(19) 
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where |a+(t)) and |a_(t)) are the displaced coherent states: 

l®±(^)) = exp i^z/a^a ± g{a + 1|q;o). (20) 

The expectation value of the position of the oscillator is given by 

{x{t)) = c{x)a_{t) + (1 - c)(x)„+(f). (21) 

The expectation value of the oscillator velocity is the quantity of interest for evaluation 
of the work flux and can be shown to be 


= 2i'sin(i/() yxc, ± ^), 


( 22 ) 


where the coupling g is also assumed to be real. The hnal expression for the work flux 
into the spin thus reads 

Wi{t) = 2g(l — 2c) [uxo + (1 — 2c)g~\ sin(z/f). (23) 

The corresponding heat fluxes are found to be, 

Qi{t) = 0, (24) 

Q 2 {'t) = 8c(l — c)( 7 ^ sin(z^f). (25) 


The heat flux into the oscillator is independent of the initial displacement of the oscillator 
and is maximum for a maximally mixed initial state of the spin (c = 1/2). The work 
flux, however, is maximized for an initially pure state (c = 1 or c = 0). The initial state 
of the spin thus serves as a control parameter that will determine the nature of energy 
exchange between the two subsystems. 


2.4- Non-commuting interaction 

This scenario is characterized by the commutators: [Ai, Hi] ^ 0, and [i? 2 ,772] A 0. From 
the effective Hamiltonian of system I, we identify Hf^{t) = 0, and Hf^{t) = Hf^{t). 
Work on system I is thus given by 

WM = (26) 

and similarly for system II: W 2 {t) = —iTr 2 {[H 2 , i7|®(f)]p2(t)}- Heat flux into system I is 
given by 

Qi{t)= -tTi{[{Hi^l2),Hi2\Ci2{t)}. (27) 

Similar expressions can be obtained for system H. The heat fluxes Qi{t) and Q 2 {t) 
are, in general, non-vanishing and unequal. As the dynamics always correlate the two 
subsystems, regardless of the choice of the initial state, Ci 2 {t) cannot be zero for all t. If 
the condition [Hi,Hi 2 ] = — [i72,77i2] is satished, however, the net heat and work fluxes 
will be zero, that is Qi{t) + Q 2 {t) = 0 and Wi{t) H 2 (t) = 0. 

The Jaynes Cummings model with the rotating wave approximation is an example 
of such an interaction: H = ^az + r'o^a-l-g^a^a-l-a-a^) where cr+ and (T_ are the raising 
and lowering operators of the spin. The condition [Hi,Hi 2 \ = —[772,77x2] translates to 
an oscillator that is on resonance with the spin, u = ujq. 


3. Generalization to open systems 


3.1. One partition becomes a Markovian bath 

In this section we demonstrate that heat flux into system I converges towards known 
results, in the limit that system II acquires the character of a macroscopic Markovian 
thermal bath. The exact equation of motion for pi{t) is given by 

pi{t) = -i[Hi,pi{t)] - iTT2[Hu, p{t)]. (28) 

Recalling that p(t) = pi(t) ® P 2 (t) + C'i 2 (t), we obtain 

pi(t) = -i[Bi,pi(t)] -iTr2[IIi2,pi(t) 0 P2(t) + Gi2(t)]. (29) 

As system II is now a thermal bath, we make the substitution P 2 (t) —>■ Pb(1) and assume 
a general system-bath interaction of the form Hib = Xi 0 Yb, such that 

Pi(t) = -^[H,+X,{YB),p,{t)]-^TTB[H,B,C^B{t)]. (30) 

For weak system-bath interactions, the expectation value (Yb) can be approximated 
by its initial thermal value. The term Xi{Yb) thus represents a bath-induced 
renormalization of the system energy. For brevity we assume (Yb) = 0. Bath operators 
are often defined to have vanishing fluctuations such that the influence of the bath can 
be entirely characterized through a two-time correlation function [23]. This is true for 
a thermal bath, provided that [i^B, Yb] ^ 0. The case of constant (Yb) renormalizes 
the system Hamiltonian and should be accounted for in the derivations that follow. 
Consistency demands that the second term in Eq. fl30|) approaches a Lindblad dissipator 
in the Markov limit: 

-tTrB[H^B,C,B{t)] A[pi(t)]. (31) 

This indicates that, even in the Markov limit, the existence of system-environment 
correlations is crucial for transfer of heat. The heat flux then takes the expected form 
in this limit, i.e. Qi{t) = Tri|ifi£i[pi(t)]}. 

Eq. (j3T]) implies that system-bath correlations are not absent in the Markov limit. 
Note that there is no contradiction between this observation and the Born approximation 
used in the derivation of Markovian master equation, which assumes separability of the 
state of the system and the bath, p pe pg(t) 0 pf,. Indeed, the Born approximation 
does not imply that there are no excitations in the bath caused by the system [23] . 
The combination of the Born and the Markov approximations provide a description on 
a coarse-grained time scale. The assumption is that environmental excitations decay 
on a time scale that is much faster than the time scale of the system dynamics, and 
thus cannot be resolved. Over this coarse-grained time scale, the state of the bath 
can be assumed to be approximately constant, and the system-bath correlations can be 
neglected. In actuality, build-up of correlations between the system and the bath leads 
to loss of phase information from the system. Care therefore must be taken in applying 
the Born approximation outside the scope of the derivation of the Markovian master 
equation. 



9 


3.2. One partition interacts with a Markovian hath 

The notion of a quantum system that serves as a time-dependent drive for a second 
quantum system, can still be applied if either, or both, partitions interact with a bath. 
In this section we present this generalization to open quantum systems. Again, note 
here the semantic difference with classical thermodynamics where such system would 
be denoted as ‘closed’ indicating that it may exchange energy (but not matter) with a 
bath. Consider an open bipartite system, where only one partition (always I) interacts 
with a Markovian bath. The Hamiltonian of the total system, including the bath, can 
be written as 



H = H^ + H2 + Hi2 + Hb + Hib, 

(32) 

where Hb 
partition I. 
form Ci[pi 

is the bath Hamiltonian and Hib describes the interaction of the bath with 
The action of the bath is described by a Markovian dissipator in the Lindblad 
], and the equations of motion for pi(t) and p 2 (t) now read 


pi{t) = -i[Hi + H^^{t),pi{t)] + Cf^[p{t)] + Ci[pi], 

(33) 


P2it) = -i[H2 + Hf-{t),p2{t)] + Ct%it)]. 

(34) 

The heat flux to system I now has an extra contribution due to the bath: 



Q.(«) = 

(35) 


02(«)=Tr2{ir'(«)£f[p(i)]}. 

(36) 


We consider the case where [H[{t),HiB\ 7 ^ 0, such that the bath can exchange energy 
with system I, and investigate the modification to the work and heat flux relations 
derived in section 2. We will not discuss the heat flux from system I to the bath, but 
note that, unless for special cases, its magnitude cannot easily be deduced from Q\{t). 

3.2.1. A fully commuting interaction: [Hi,Hif\ = 0, i = {1,2}. Unlike the 
corresponding closed system scenario (Sec. 2 . 2 ), the work and heat fluxes on the 
subsystems are unequal prior to equilibration and are given by 

Wi(t) = (H 2 (t))(Ai(f)) = 0, (37) 

W2{t) = (H2(t))(ii(t)) = (H2(t))Tri{Ai£i[pi(t)]}. (38) 

Interestingly, the bath-induced relaxation of the operator (Ai(t)) is manifested as a work 
flux on system II until (Ai(f)) = 0 at thermal equilibrium. As for heat we obtain 

Qi(f) = Tri{i7((f)£i[pi(f)]}, (39) 

W) = 0. (40) 

The heat flux Qi(f), is maintained until system I is thermalized, at which point 
Qiit) = Q 2 {t) = 0. The dispersive spin-oscillator Hamiltonian with the spin exposed to 
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a heat bath provides an example of this scenario: H = pa'a + ujQaz+ga^aaz + HB + Hi-q, 
where Hb = and Hib = ^ J2k^k{bk + &!)• The interaction of the spin 

with the bath will result in (i) a net transfer of heat to the spin and (ii) coherent energy 
transfer from the spin to the oscillator. For instance, the oscillator can represent the 
mechanical mode of a resonator, the spin a two-level atom, and the bath the background 
radiation. The relaxation of the atom is subsequently manifested as a net work flux on 
the resonator. 

3.2.2. A partially commuting interaction. Given that in our analysis the bath always 
interacts with system I, we consider two different scenarios of a partially commuting 
interaction to cover all possibilities. 

Case 1. [Hi, = 0, [H 2 , ^^ 12 ] 7 ^ 0. Work fluxes in this case are given by 

Wi{t) = -W2{t) = {Ai){B2). (41) 

The condition Wi{t) + W 2 it) = 0 is thus still satisfied (cf. Eq. ffTS]) ) despite the addition 
of the thermal bath. The heat fluxes are modihed as follows, 

Qi(i) = Tli{ff;(i)Abi(i)|}, (42) 

4(4)= -iTr{|(l, (43) 

In the corresponding closed scenarios, we observed that heat flux is solely into system II. 
This observation is no longer valid for the open system prior to thermal equilibrium. 

An example of such an interaction is provided by the displaced spin-oscillator 
Hamiltonian with the spin exposed to a bosonic bath; H = -|- gczia^ -l- 

a) + Hb + Hib, where Hib = (Tx<S)J2k ^k{bk + b\). 

Case 2. [Hi,Hi 2 \ 7 ^ 0, [H 2 ,Hi 2 \ = 0. The work flux takes the same form as the 
previous case, but the heat flux is now exclusively into system I: 

Qi{t) = -*Tr{[(li ® i/'(t)),i/i 2 ]C'i 2 (t)} +Tri{i7((t)A[pi(t)]}, (44) 

Q2{t) = 0 (45) 

An example of such a system is provided by the spin-oscillator model where now only the 
oscillator interacts with a bosonic bath, that is, H = ‘^az + i'a^a+gaz{a^+a)+HB+HiB, 
for Hib = {a + a^) 0 Y.k ^k{bk + bl). 

3.2.3. A non-commuting interaction. For a fully non-commuting interaction, work 

. nn 

flux associated to each subsystem is given by Wj{t) = —iTij{[Hj,Hj^{t)]pj{t)}. The 
heat fluxes take the form 

Qiit) = -zTt{[{Hi 0 l 2 )),ifi 2 ]Gi 2 (t)} +Tri|ih((f)A[pi(f)]|, (46) 

Q2{t)= -tTT{[{li^H2),Hi2]Ci2{t)}. (47) 

Similar to the hnite system case (Sec. 2.4.), if [Hi,Hi 2 ] = —[H 2 ,Hi 2 ], then Wi(t) -\- 
W{t) = 0. At thermal equilibrium, where the heat flux due to the bath vanishes, this 
condition also implies that Qi(t) -|- Q(t) = 0. 
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The Jaynes Cummings model is an example of such an interaction: H = ^cTz + 
ua^a + 5 '(cr+a + cr_al) + Hb + Hib- The bosonic bath can interact with either the spin 
[e.g.: Hib = crx<S)J2k ^k{bk + bl)] or the oscillator [e.g.: Hib = (a + al)(8)Xlfc Afc(6|, + 6|,)]. 

3.3. Each subsystem interacts with a bath 

The present treatment can readily be extended to the case where both sub-systems 
interact with a thermal bath. We will only consider the case of the fully-commuting 
interaction for brevity, but the theory can straightforwardly be applied to all cases. 
Consider a bi-partite system where each partition interacts with its own Markovian 
bath. The total Hamiltonian is given by 

H = H, + H2 + Hu + Hb, + Hb, + H^b + H 2 B, (48) 

where Hb, and Hb, are the bath Hamiltonians, while Hib and H 2 B specify the system- 
bath interactions. The action of the two baths is described by the Markovian dissipators 
£i[pi] and C 2 [p 2 ]- The equations of motion for the partitions now take the form 

Pj (t) = -i [Hj + Hf^{t) , pj (t)] + ® [p(t)] + Cj [pj]. (49) 

For the fully commuting interaction, work fluxes are given, Wi = {B 2 {t)){Ai{t)), W 2 = 
{B 2 {t)){Ai{t)), while the heat flux takes the form Qj{t) = Trj|Trj(f)£j[pj(t)]}. In the 
steady state, energy conservation and vanishing work fluxes mean that all energy fluxes 
are zero. For the fully non-commuting interaction, however, it can be shown that a 
hnite heat and work flux can be maintained in the steady-state. Even with a fully 
commuting interaction, more elaborate designs consisting of multiple baths of different 
temperatures can be utilized to engineer thermodynamic devices capable of maintaining 
a constant heat flux in the steady state. For instance, there are prototype thermoelectric 
devices consisting of two interacting quantum dots coupled to three thermal reservoirs 
|24] . The coupling between the quantum dots is a Coulombic repulsion and therefore 
a fully-commuting interaction. It has been demonstrated both experimentally pa and 
theoretically [ 23 ] that this design aids the extraction of energy from a temperature 
difference while avoiding any direct connection between the hot and cold reservoirs 
along which heat might flow. 

4. Entropy change and the second law 

In Ref. [1^ it is argued that entropy and heat flux in a hnite quantum system are 
related via the expression 

dSi = /3*dQi, (50) 

where dSi is the change in the von Neumann entropy of system I and f3* is dehned to 
be an effective inverse temperature. This dehnition is inspired by analogy with classical 
thermodynamics and was used to demonstrate agreement with Clausius relation for 
systems at thermal equilibrium. There is, however, a problem with this identihcation as 
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counterexamples to Eq. fl5U]) may be found where the heat flux is vanishing {dQi = 0), 
but the entropy flux is hnite {dSi 7 ^ 0). To illustrate this point, consider the 
scenario discussed in section 2.2 where both subsystem Hamiltonians commute with 
the interaction. In this case, there is no heat flux into either of the partitions, yet the 
local entropy of system I can still change, provided that it has an initial coherence and 
that partition II does not start in an energy eigenstate. In appendix A we show that 
if these two conditions are satisfied, the evolution of partition I is non-unitary, and its 
local entropy is oscillatory, despite the lack of a heat flux, and in contradiction with Eq. 
dSOD. 

Clearly, a different line of thought ought to be followed to reconcile the energy flux 
formalism with the second law. For an initial separable state p(0) = pi(0) $§ P 2 ( 0 ), 
evolving under the Hamiltonian H = Hi + H 2 + i?i 2 , Esposito and et ah [26] have 
demonstrated that the change in the von Neuman entropy of system I, between t = 0 
and t = T, can be expressed as 


A^i(r)=^i(r)-^i(0) 


= D 


P(^)l|Pl(^)P2(0) 


+ Tr<^ [p2(r) -p2(0)]lnp2(0) \ , 


(51) 


where Zl[pi||p 2 ] = Tr{pilnpi} — Tr{pilnp 2 } is the relative entropy. This expression 
is derived from the conservation of the total entropy with the sole assumption of the 
separability of the initial states. The first term in Eq. (j5Ti) is the irreversible contribution 
to the entropy change and is known as the entropy production. The second term is the 
reversible contribution to the entropy change. It quantifies the exchange of heat and 
results from the reduced dynamical evolution of system H. The entropy change can 
therefore be written in the standard thermodynamic form AS'i(r) = + AS*^® with 


A^i’'(r) = D 


/>(t)IIpi(t)p2(0) 


AS;*(t) = Tri [p 2 (t) - P2(0)]liip2(0) 


(52) 

(63) 


Since the relative entropy is positive, Eq. fl52l) indicates the positivity of the entropy 
production, i.e.: AS'J'’(r) > 0 , thus fulhlling the second law of thermodynamics [ 2 ^ [27j . 

Equations (j52l) and (153|) can also be used to investigate any possible links between 
heat and entropy. Clearly, Eq. fl53l) cannot be identified with a heat flux. In classical 
thermodynamics the association between heat and entropy holds for reversible entropy 
exchanges and is defined for a thermal initial state. Eq. fl52l) shows that this associa¬ 
tion is not correct for the general quantum setting. The irreversible entropy is due to 
correlations between the two partitions that may arise independently of a heat flux. We 
therefore conclude that the direct association between heat flux and entropy change, as 
suggested by Eq. fl50l) . does not hold in general. 


As an example, equations fl52|) and fl53l) can be used to analyze the entropy flux for 
the commuting interactions (sec. 2.2.). From Eq. (j52i) . we conclude that AS'{®(r) = 0, 
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Figure 2. (a) Heat and work contributions to the total energy change AUi{t) of 
a two level system as it interacts with a single harmonic mode, (b) Reversible and 
irreversible contributions to the change in the von Neumann entropy of the two level 
system. Parameters are T = 300 K, ujq = 150 cm“^, H = 50 cm“^, ly = 180 cm“^, 
and g = 50 cm“^. 

provided that system II is initially in a statistical mixture of its energy eigenstates. The 
heat flux is also zero for commuting interactions. Yet, if system I possesses some initial 
coherence, the dynamics would correlate the two partitions and the entropy production 
AS'l^(r) would exhibit oscillations. 

To explore further possible relations between AS'[®(r) and AQi(r), we consider the 
non-commuting case where a two-level system interacts with a single harmonic mode. 
This spin-boson type Hamiltonian is given hy H = + Va^) + z/a^a -|- g<Jz{ci + al), 

with bias Wq, tunnelling amplitude V, mode frequency n and system-mode coupling 
strength g. The spin is initially assumed to be in the excited state while the oscillator 
is initially in a thermal state at a temperature T. For these initial conditions the 
reversible entropy change quantifies the energy change of system II and becomes 
= —I32{{H2)t — {H 2 )o) with (32 = 1/T. Notice that this energy change for 
system II is not trivially related to the heat flux into system I. Fig. |2^ displays the 
heat and work contributions to energy changes of the spin (parameters provided in 
the caption of Fig. [2]), indicating there is predominantly a negative heat flux into 
the spin i.e. the system is “cooled”. The reversible and irreversible contributions to 
entropy change of the spin are displayed in Fig. [2b. ASI^{t) is negative indicating 
that as the spin and oscillator become correlated, the latter evolves towards a non- 
thermal state with an average energy larger than its initial energy. Although AS'J®(r) 
resembles some features of AQi(r), such as its negativity and the position of some 
of the extrema, a proportionality relation of the form of Eq. (8) in Ref. [26] does not 
follow straightforwardly. We conclude this section by noting that although entropy 
production is positive its derivative can be negative. This is a direct consequence of the 
non-Markovian nature of the reduced dynamics of system I, and positivity of the rate 
is recovered in the Markov limit [26] . 

In the section, we have demonstrated that the energy flux formalism satishes the 
second law of thermodynamics. By partitioning entropy changes as proposed in Ref. 
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[25], the positivity of entropy production within the heat flux formalism is demonstrated. 
The direct relation between heat and entropy proposed in Ref. m is incorrect in general. 

5. Conclusions 

In this article we have considered internal exchange of energy in bipartite quantum 
systems and have investigated the consequences of the partitioning introduced in Ref. 
[HI for different commutations of the interaction and the bare Hamiltonians. A 
generalization of the formalism to Markovian open quantum systems was also presented 
and the consistency of the formalism with the second law of thermodynamics was 
discussed. 

From a classical standpoint exchange of heat and work have distinct physical origins; 
work is exchanged due to the action of an external drive. Heat, on the other hand, is 
exchanged with a thermal bath. The present formalism connects these two scenarios 
and states how each situation arises in a particular limit of the general bipartite system. 

Optomechanical implementations provide the ideal platform for the investigation 
of energy transfer between quantum mechanical systems j 28 | and the prototype 
spin-oscillator models can readily be implemented in such a setting. Nano-scale 
thermoelectric devices comprised of quantum dots interacting with thermal reservoirs are 
alternative systems in which the interplay of quantum mechanical and thermodynamic 
principles can be investigated pa EH]. There are, however, still issues associated with 
the measurement of the aforementioned fluxes and the practical applicability of the 
energy partitioning on the entropy fluxes requires further investigation. 
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8. Appendix A 

In this appendix we demonstrate that the reduced dynamics under a fully commuting 
interaction can be non-unitary. We consider the Hamiltonian H = Hi + H 2 + H 12 where 
[Hi, Hi 2 \ = [H 2 , Hi 2 \ = 0. The interaction Hamiltonian can be written as H = Ai ^ B 2 
where Ai acts on partition I and B 2 on partition H. Hi and Ai can be expanded in the 
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energy eigenbasis of system I. That is, 

Hi = '^ei\ei){ei\, (54) 

i 

Ai = '^ai\ei){€i\, (55) 

i 

where {|ej)} are the eigenstates of system I with eigenenergies {ej}, and {oj} are a set 
of coefficients. Similarly, we expand H 2 and B 2 in the eigenbasis of H 2 

H 2 = ^Ek\Ek){Ek\, (56) 

i 

B2 = Y,bk\Ek){E,\, (57) 

k 

where {\Ek)} are the eigenstates of system II with eigenenergies {E^} and {b^} 
are a set of coefficients. The reduced dynamics of partition I is given by pi{t) = 
Tr 2 { exp(—ii/t)p(0) exp (iTTt)}. We expand the Hamiltonians and arrive at the 
following hnal expression 

pi{t) = (g) (E,|p2(0)|^fc). (58) 

ijk 

Unless partition II starts in an energy eigenstate, the dynamics of coherences {ei\pi(t)\ej) 
are non-unitary. We therefore conclude that pi(t) exhibits non-unitary dynamics 
provided that system II is not in an energy eigenstate at t = 0, and system I 
possesses some initial coherence. Non-unitary dynamics give rise to a time-dependent 
von Neumann entropy for system I. 


